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Abstract 

Exact global solutions of five-dimensional cosmological models compactified on a S\/ Z2 orbifold 
with two 3-branes are presented, and evolution of a simple model is studied. It is found that on 
all 4D spacetime hyper surf aces, except a singular one, the expanding universe was started not 
from a big bang but from a big bounce, and before the bounce the universe was in a deflationary 
contracting phase. It is also found that whether the energy density on the second brane is positive, 
zero, or negative depends on the size of the fifth dimension. 
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It was proposed [1-3] that our universe is a 3-brane embedded in a higher dimensional 
space. While gravity can freely propagate in all dimensions, the standard matter particles 
and forces are confined to the 3-brane. Binetruy, Deffayet and Langlois (BDL) have consid- 
ered a 5D cosmological model and derived the Friedmann equations on the branes [4]. This 
model has received extensive studies and most of these studies were focused on the branes 
only [5] . Be aware that it is technically difficult to move from the branes and construct exact 
solutions of the field equations in the bulk, whereas numerous exact solutions are known in 
5D Kaluza-Klein theories. This suggests that we can take a solution of the Kaluza-Klein 
theory and use the Z2 reflection symmetry and Israel's jump conditions to obtain a global 
brane model. In what follows, we will do this for a rich class of 5D cosmological solutions 
found originally by Liu and Mashhoon [6] and restudied recently by Liu and Wesson [7,8]. 
We will then construct a class of global brane solutions which is of the simplest BDL'type, 
i.e., the matter on the two branes is a perfect fluid, the bulk is empty, and no cosmological 
constants are added in the bulk and on the branes. 

The 5D metric for the solutions is 

dS 2 = B 2 dt 2 -A 2 ( + r 2 dn 2 ) - dy 2 . (1) 
\ 1 — kv J 

where B = B{txf) and A = A(t t y) are two scale factors, k(= ±1 or 0) is the 3D curvature 
index, and dVl 2 = (d6 2 + sin 2 6d(f) 2 ). A class of general solutions of the 5D vacuum equations 

R AB = 0, A,B — 0123; 5 (2) 

was given [7] by 

,2 / 2 ,\ 2 n v 2 + K ^ IdA A 

where an overdot denotes partial derivative with respect to t, fi — fj,(t) and v = u(t) are 
two arbitrary functions, and K is a 5D curvature constant related to the square of the 
Riemann-Christoffel tensor (Kretschmann scalar) via 

RabcdR ABCD = ^ 8 ■ (4) 

To confirm solutions (jSJ), one can go to see the detailed derivation of the solutions in Ref. [6] 
(in which different notations were used); or one can substitute (JSJ) into Rab directly (which 
is complicated); or one can use computer programs such as MAPLE or GRTensor. 



Because solutions (J3J) satisfy the 5D vacuum equations (0 ), we can use them as the 
bulk solutions of the BDL-type brane models. Note that one can change y to —y without 
violating the validity of (0) as exact solutions of (j2J)- Therefore, to obtain brane models we 
use the Z 2 reflection symmetry on (jHJ), i.e., we let 

,2 f 2 i\ 2 n i i v 2 + K IdA A 

Then we take the corresponding 5D Einstein equations being 
Gab = ^(5)^5 

+% - y 2 ) diag (p 2 , ~P2, ~P2, -P2, 0) . (6) 

Here the first brane is at y = y\ = and the second is at y = y 2 > 0. In the bulk we have 
Tab = and then Gab = 0, so equations (jOJ) are satisfied by (j^J). On the branes we have to 
solve the equations (0) as follows. 

We use the 5D metric (JTJ) to calculate the 5D Einstein tensor, which, substituted in ©, 
gives 

g »- 3 (I + 4)- 3b2 (t+S 

= Af 5) B 2 [ Pl %) + p 2 5(y - sj)] , (7) 
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According to Israel's jump conditions, the two scale factors A and B are required to be con- 
tinuous across the two branes. Their first derivatives with respect to y can be discontinuous 
across the branes, and then their second derivatives give a Dirac delta function. So we have 
to calculate A' and B' across the two branes. By differentiating A 2 in (JSJ) with respect to y, 
we obtain 

AA'=(^ + k)y-u?M. (11) 



dy 



Therefore we get 



where A\ = A(t,y — y\ = 0) and A 2 = A(t,y = y 2 ). Then, using (JT2J) and the second 
equation in (JHJ), i.e., B = A/fi, we get 



H&t \Ai J fidt \A 1 

, s 1 d f u 2 + k v \ _.. , I d { u 2 + k , > 

= ( V s2 - • = ^ ( V* - • (13) 

So the jumps of A' and 5' across the two branes are 

r „n 2l/ r ... f U 2 + k V 

|B| ' = -^UJ- l B l> = -JIa(— »-aJ • (14) 

where [A'^ = v4'(0 + ) — A'(0~) and so on. Substituting (fT4^1 in the field equations (Ej) and 
fllOjl . we obtain 

(5)Pl = ~A [ ]l = A 2 ' 

^4^4^=4=^4, us) 

and 



3 6 f n 2 + k v 

1 r „„ 2 



-2/2 - -r - -r —^2/2 - -r • (16) 
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Thus we have derived the energy densities and pressures on the two branes. Meanwhile, the 
conservation law Tj* B = gives 

Pi + 3(Pi+Pi)j; = 0, i = l,2. (17) 

This relation can also be verified directly by substituting ()15|). ( HSj) and (j5J) into (fT7)l as 
expected. 

From the 5D metric (JTJ we see that on a given y = constant 4D hypersurface the proper 
time can be defined as dr = B(t,y)dt. So the Hubble and deceleration parameters can be 
defined as 

^>4M> ■ft'>--£!®/®'--2 M 

where we have used the relation B = Aj\i. Meanwhile, the first equation in (jSJ for the 
branes can be written as 

(fi 2 + k)A 2 = [( fJ 2 + k)y i -u] 2 + K , i = l,2. (19) 

Then, using (fTKjl . (fTCl) and <jTB)l in <fTTTjl - we obtain 

# + 4 = ^M+£> * = (2°) 



^ ( i_ 2? .) + _ = _J| p .( p . + 2 p.)-_, < = 1,2 . (21) 



These two equations are the induced Friedmann equations on the branes. 

Thus we obtain a complete set of global exact solutions given in (JHJ) and (|15j) - (j21j) which 
contains two arbitrary functions fi(t) and v(t). From the relation B = Aj \i and metric 
(Q) we see that the form of Bdt is invariant under an arbitrary coordinate transformation 
t — > t(t) .. This freedom can be used to fix one of the two functions /x(t) and u(t). Another 
freedom corresponds, as is in the standard general relativity, to the unspecified equation of 
state of matter. For a given equation of state p = p(p), the two equations in (fTK|) give a 
constraint that can be used to determine the function, say, u(t). So, generally speaking, 
if the matter content on the first brane is known, then fi(t) and v(t) can be fixed. Then 
the whole solutions can be fixed too. Then, by ()16j) . if the size of the fifth dimension 
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is also known, we will know the matter content on the second brane, and this is of great 
importance. 

Using above procedure, we can, in principle, construct particular models such as the 
matter and radiation dominated models as in the standard FRW cosmology. However, 
preliminary studies show that they are not mathematically simple. So we leave this study 
in the future. Here, in this letter, we will use another way to obtain an explicit universe 
model, that is, we will choose fi(t) and v(t) firstly and then to study the evolution and 
matter properties. Now we let 



k = 0, K = ^, u(t) = j- t , 



n(t) = {2Lty 1/2 



(22) 



where L is a constant with a dimension of length and t b (we assume t b > 0) is a critical 
constant with a dimension of time. In this way the solutions © become 
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So on the first brane we have 
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On the second brane we have 
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4t 2 _ (2/2 _ 2tb y 4t 2 + (2/ 2 - 2t b y 

Equations (|2*S|) - (|2*3|) constitute a quite simple two-brane model. 

From (|23|) we see that, in a y = constant hypersurface, the scale factors B m 1 and 



A w ^J2t/L for 2t ^> — 2tb\. So it approaches the standard radiation-dominated model 
at late times of the universe. The global evolution of the scale factor A(t, y) with t b = 1 
and L = 1 is plotted in Fig.l. From this figure we see that there is a singular surface 
\y\ — 2t b in which A —>■ as t —>■ , showing a big bang singularity as is in the standard 



FIG. 1: Global evolution of the scale factor A(t, y) = y2t + — 2) 2 / (2t). The first brane (our 
universe) is at y = 0; while the second brane (the hidden universe) could be placed at y = 2, 
< y < 2, or y > 2, giving three types of brane cosmological models, respectively. 

FRW models. However, in all other hypersurfaces, the scale factor A reaches a non-zero 
minimum y2t/Z at 2t = \\y\ — 2t&| (where B = 0); and at both sides of this minimum, 
A tends to infinity. So this minimum can naturally be explained as a big bounce [7]. 
Generally, the solution has two kinds of singularities, corresponding to A = and B = 
respectively as discussed in Ref. [7]. A = is, by (j3J, an intrinsic singularity of the 5D 
manifold, whereas B = is just a coordinate singularity similar to the Schwarzschild event 
horizon. So here the big bounce singularity belongs to the second kind. 

The first brane is at y — 0. By (|24p we see that at the bounce point t = tb we have Ai 
tends to its minimum, Bi tends to zero, p l is finite, and p± tends to infinity. Physically, we 
can say that it is this infinitely large pressure pi that caused the big bounce. 

From 125(1 we see that the evolution and the matter properties on the second brane depend 
on the size 2/2 of the fifth dimension (see Fig. 1 also). Thus we have three types of brane 
models as shown below. 

Type I. The second brane is at y 2 = 2t b . Then (J23j) gives A 2 = ^/2t/L, B 2 = 1, and p 2 = 
p 2 = 0. This implies that the second brane is empty of matter. So Type I is actually 
a one-brane model for which there was a big bounce on our side of the model and 
a big bang on the "hidden" side. 



Type II. The second brane is at < 1/2 < 2£&. Then (|25|) shows that there were bounces 
on both branes. On the second brane we have p 2 < 0. So Type II is a two-brane 
model for which the energy density on the second brane is negative. 

Type III. The second brane is at y 2 > 2%. This is also a two-brane model for which the 
energy density on the second brane is positive. Specifically, if y 2 = 4t&, then A 2 = A\, 
B2 — Bi, p 2 = Pi, and P2 — Pi, giving a completely symmetric two-brane model. 

Now let us consider the solution (J24j) for which when the coordinate time t varies from 
to tb and then to +00, the spatial scale factor Ai(t) contracts from +00 to a minimum 
Anin = 2\/tb/L > and then expands to +00 again. We also find that as t — > 0, A\ — > 
2t\j[Lt\ B\ — > (t b /t) 2 and the 4D metric on the first brane tends to 




By a coordinate transformation 

t = 2L-H 2 b e T/tb , (27) 

we get 

ds\ -> dr 2 - (dr 2 + r 2 dQ 2 ) (as r -> -00) , (28) 

which represents a deflationary de Sitter cosmological model. Note that by (|2*7j) t — 
corresponds to r = —00. Therefore we conclude that according to the proper time r, the 
universe on the first brane has been existed forever and contracts from a de Sitter space to 
a non-zero minimum, at which the pressure pi reaches to infinity and causes a big bounce. 
After then, the universe expands. 

In conclusion, we have derived a class of exact global solutions of five-dimensional cos- 
mological models with two 3-branes, for which matter on the branes is of the form of a 
perfect fluid, the bulk is empty, and no cosmological constants were introduced. This solu- 
tions contain two arbitrary constants k and K, corresponding to the 3D and 5D curvatures 
respectively, and two arbitrary functions p(t) and v(t), corresponding to the two freedoms: 
the arbitrary coordinate transformation t — > t(t) and the unspecified equation of state of 
matter. By choosing p(t) and v(f) properly a special model is discussed in detail. It is 
found that in the late times of the universe this model evolves at the same rate as is in the 
standard radiation-dominated FRW model. And there is a singular 4D hypersurface y = 
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on which there was a big bang. On all other hypersurfaces there were no big bang but big 
bounces. Before the bounces, the universe was in a contracting phase. It is also found that 
whether the energy density on the second brane is positive, zero, or negative depend on 
the size of the fifth dimension. Thus three types of models are obtained for which type I is 
actually a one-brane model while types II and III are two-brane models with negative and 
positive energy densities, respectively. 
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